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Abstract 



We compute the baryon asymmetry generated at the electroweak phase tran- 
sition by the Higgs scalar sector of the minimal supersymmetric standard 
model. Because of large enhancement effects from low momentum modes, 
Higgs particles may be responsible for the observed baryon asymmetry even 
though CP-violation in the Higgs sector only appears at the one-loop level. 
We also discuss the approximations made in the analysis and suggest possible 
improvements. 
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1. The Standard Model (SM) fulfills all the requirements for a successful generation of the 
baryon number at the electroweak scale due to the presence of baryon number violating 
processes which also impose severe constraints on models where the baryon asymmetry is 
created at very high energy scales 0. Unfortunately, the electroweak phase transition is at 
best weakly first order in the SM meaning that the baryon asymmetry generated during 
the transition would be subsequently erased by unsuppressed sphaleron transitions in the 
broken phase. Therefore, if one is willing to pursue the idea of electroweak baryogenesis, 
some new physics at the weak scale must be called for. The most promising and well- 
motivated candidate seems to be supersymmetry (SUSY). Electroweak baryogenesis in the 
framework of the Minimal Supersymmetric Standard Model (MSSM) has attracted much 
attention in the past years with particular emphasis on the strength of the phase transition 
10-0 and the mechanism of baryon number generation [0-0. 

Recent analytical |jlO[ and lattice computations |ll| have revealed that the phase tran- 
sition can be sufficiently strongly first order if the lightest stop is not much heavier than 
the top quark, the ratio of the vacuum expectation values of the two neutral Higgses tan j3 
is smaller than ~ 4 and the lightest Higgs is lighter than about 85 GeV. Moreover, the 
MSSM has the bonus that additional sources of CP-violation may be present besides the 
CKM matrix phase. These new phases are essential for the generation of the baryon number 
since large CP-violating sources may be locally induced by the passage of the bubble wall 
separating the broken from the unbroken phase. Baryogenesis is fueled when transport prop- 
erties allow the CP-violating charges to efficiently diffuse in front of the advancing bubble 
wall where anomalous electroweak baryon violating processes are not suppressed. The new 
phases appear in the soft supersymmetry breaking parameters associated to the stop mixing 
angle and to the gaugino and neutralino mass matrices; large values of the stop mixing 
angle are, however, strongly restricted in order to preserve a sufficiently strong first order 
electroweak phase transition. Therefore, an acceptable baryon asymmetry from the stop 
sector may only be generated through a delicate balance between the values of the different 
soft supersymmetry breaking parameters contributing to the stop mixing parameter, and 
their associated CP-violating phases [||. On the other hand, charginos and neutralinos may 
be responsible for the observed baryon asymmetry if the phase of the parameter /i is larger 
than about 0.1 [|,||. 

Since in the case in which the bubble wall is thick, the Higgs background carries a 
very low momentum (of order of the inverse of the bubble wall width L^), the production 
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of the baryon asymmetry is enhanced if the degrees of freedom in the gaugino/neutrahno 
sectors are nearly degenerate in mass ||^: values of the phase of fi lower than 0.1 are only 
consistent with the observed baryon asymmetry for values of of order of the gaugino 
mass parameters. This is due to a large enhancement of the computed baryon asymmetry 
for these values of the parameters. 

The fact that CP-violating sources are most easily built up by the transmission of low 
momentum particles over a distance is an indication that particles with masses 

much smaller than the temperature T may be relevant in the process of quantum interfer- 
ence leading to CP-violating currents in the bubble wall. It is therefore fair to guess that 
the nearly massless neutral and charged modes present in the scalar Higgs sector during 
the electoweak phase transition may play a significant role in supersymmetric electroweak 
baryogenesis. This is analogous to what happens in the SM where the barrier in the effective 
potential separating the two minima at finite temperature is originated by infrared effects. 

CP-violating Higgs currents are generated by the interactions with the bubble wall only if 
some CP-violation is present in the scalar potential. This does not happen at the tree-level 
because the form of the scalar potential is dictated by supersymmetry and CP-violating 
operators are forbidden. However, once supersymmetry is broken and if CP-violation is 
present in the stop and in the -ino sector, CP-violating Higgs operators are induced at 
finite temperature by one-loop diagrams involving stops, charginos and neutralinos [IT2 



CP-violating scatterings of the Higgs modes with the adavancing bubble wall are therefore 
expected to give rise to some baryon asymmetry. 

The goal of this Letter is to compute the final baryon asymmetry generated by the Higgs 
scalar sector. Because of the large enhancement from the infrared region, Higgs particles 
may be responsible for the observed baryon asymmetry even though the CP- violation in the 
Higgs sector only appears as a one-loop effect. The details of the computation are contained 
in the next section, while section 3 is devoted to our conclusions and comments. 

2. To compute the CP- violating Higgs current we make use of the method proposed 



in Ref. |T3[ and recently adopted by Carena et al. to compute the stop and the -ino 
contribution to the baryon asymmetry. This method is entirely based on a nonequilibrium 
quantum field theory diagrammatic approach. It may be applied for all wall shapes and 
sizes of the bubble wall and it naturally incorporates the effects of the incoherent nature 
of plasma physics on CP-violating observables. What we need to compute is the temporal 
evolution of a classical order parameter, namely the CP- violating Higgs current, with definite 
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initial conditions. In this respect, the ordinary equihbrium quantum field theory at finite 
temperature may not be applied, since it mainly deals with transition amplitudes in particle 
reactions. Therefore, the closed-time path formalism is used, which is a powerful Green's 
function formulation for describing nonequilibrium phenomena in field theory . 



Following I^J^ , we are interested in the generation of some charge which is approximately 
conserved in the symmetric phase, so that it can efficiently diffuse in front of the bubble where 
baryon number violation is fast, and non-orthogonal to baryon number, so that the genera- 
tion of a non-zero baryon charge is energetically favoured. A charge with these characteristics 
is the Higgs charge density (J^) = (J^^ - J^,), where J^^j^z) = i [Hjd^H, - d^HjHi) with 
i = 1,2 and Hi and if 2 denote the two Higgs doublets. The CP- violating sources 'yniz) (per 
unit volume and unit time) of the Higgs charge density Jfj associated to the current Jfj and 
accumulated by the moving wall at a point of the plasma can then be constructed from 
Jh{z) using the relation 'yniz) = do{Jfj{z)). This is appropriate to describe the damping 
effects originated by the plasma interactions, but does not incorporate any relaxation time 
scale arising when diffusion and particle changing interactions are included. However, one 
can leave aside diffusion and particle changing interactions and account for them indepen- 
dently in the rate equations. This is a good approximation (at least for small bubble wall 
velocity) since, for instance, the typical diffusion time ~ -D/f^, where D is the typical 
diffusion constant, is much larger than any other time scales 

The CP-violating part of the MSSM Higgs scalar potential reads 

Vcp = A5 {HiH2f + Ae \Hif H^H^ + A7 li/al' HiH2 + h.c. (1) 

As mentioned in the introduction, the coefficients As^ej are zero at the tree-level, but are 
nonvanishing at the one-loop level at finite temperature when interactions of the Higgs fields 
with charginos and neutralinosQ are taken into account |12|. Since these coefficients appear 



only once supersymmetry is broken, their values are proportional to the gaugino masses 
M12, and to the supersymmetric parameter /i. The typical value for As,? is about 10~^, 
while A5 is at most 10"'* |12|. What is relevant for us is that the A's are complex if the 
parameter carries a phase 0^, A5 = lAsje*^'^'', Aej = lAejje*'^''. It is the presence of this 
phase that gives origin to the baryon asymmetry. 



^The stop contribution is suppressed when the left-handed stop mass is much larger than the 



temperature. This condition is necessary for the phase transition to be strong enough |1C]. 
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The interesting dynamics for baryogenesis takes place in a region close to or inside the 
bubble wall and we approximate it with an infinite plane traveling at a constant speed Vw 
along the z-axis. 

The computation of the CP- violating Higgs current is done by expanding (J^) around 
the symmetric phase defined by {H^) = vi{T) = and (-^2) = "^^iT) = 0. This is equivalent 
to a diagrammatic approach where (Jh) is expressed as a sum of different Feynman diagrams 
with increasing powers of external Higgs insertions. The external Higgs background has to 
be identified with the bubble wall configuration and the expansion is only justified in a 
region of the parameter space for which the mean free time r is smaller than the scale of 



variation of the masses, i.e. the wall thickness L^. The reader is referred to refs. |13 



for more details. Since the coefficients of the expansion are computed in the symmetric 
phase, we must deal with the resummation of the propagators of the Higgs fields in order 



to deal with infrared divergencies [jT5|. Indeed, to account for the interactions with the 
surrounding particles of the thermal bath, particles must be substituted by quasiparticles 
and dressed propagators are to be adopted. Self-energy corrections at one- or two-loops to 
the propagator modify the dispersion relations and introduce a nontrivial damping rate F 
due to the imaginary contributions to the self-energy. In what follows we shall adopt dressed 
propagators to compute the thermal averages of the composite Higgs operator, which allows 
us us to naturally and self-consistently include the effects of the incoherent nature of plasma 
physics. 

In the unbroken phase, the Higgs spectrum contains two complex electrically neutral 
fields and two charged ones. At the tree-level, the squared masses of one of the neutral 
states and one of the charged ones are negative, since the origin of the field space becomes a 
minimum of the effective potential only after inclusion of the finite temperature corrections. 
The resummation can be achieved by considering the propagators for the eigenstates of the 
thermal mass matrix, which has positive eigenvalues given by 



2 m^(T) + ml{T) T ^ {ml{T) - ml{T)Y + 4 mjiT) 
mh,H = ^ , (2) 

where the mf(T) are the thermal corrected mass parameters of the effective potential, 
ml{T) ~ m2 + (3/8)^27-2^ ml{T) mj + {l/2)h'^^T'^, while ml{T) = m\ (+ logarithmic 
corrections in T). Here the T = coefficients of the operator |-f^2p and 

{H1H2) respectively and g and ht are the SU{2) ^ gauge coupling and the top Yukawa cou- 
pling. Correspondingly, the neutral complex eigenstates are given by 



5 



h = ceHl* + seH2. 
H = -se H',* + ce H- 



0* I ^„ u (^) 



where ce = cos 6, se = sin 6 and 6 is the critical angle identified by the fiat direction of the 
effective potential around f i = f 2 = at the critical temperature. Analogous formulae hold 
for the charged eigenstates. 

Following refs. |l3|.p[1, it is not difficult to show that {Jh{^)) g^ts contributions from 
several one-loop Feynman diagrams which are obtained by assigning in all the possible ways 
the space-time points z and x (which is the point where the external Higgs configurations 
are attached to) on the positive or negative time branches typical of the nonequilibrium 
approach [|I^. Disregarding the T = piece, the result is (in the plasma frame) 

( J^(z)) = 8 cos^ 26 lm[fXz)] I{(3, m^, F), 

f°° , f d?k e"^'"" sin^WfcM sin/3F 
/(Am^,r)=y^ du J [cosh/5^,- COS /3F]' 

f{z) = {2(A6c2 + Xrsl)vi{z)v2{z) + 2ceSe [KvI{z) + hvl{z)\ } . (4) 

In the expression above we have only included the dominant contribution coming from the 
diagrams where only nearly Higgs massless modes at the critical temperature propagate; 
uj]: = k'^ + ml, where rrih denote the mass of the degrees of freedom propagating in the loop 
whose decay rate in the plasma is F; /3 = T~^. It is easy to show that the contribution pro- 
portional to A5 is suppressed because IA5I ^ I -^6,7 1 and because heavy states must propagate 
in the relative one-loop diagrams. The integral over the " time " variable u makes evident 
the causality inherent to the nonequilibrium approach: only the processes taking places at 
times smaller than may give rise to a nonvanishing current. 

Performing the integrals in I{(3, rrih, F) and taking the limit F <^ T, we obtain the final 
expression 

(J°(;.)) ^-cos'2eim[f{z)] (5) 

Eqs. (^ and (|^) warrant some comments. First, we notice that the momentum integration 
is infrared dominated: quasiparticles with long wavelengths and momentum perpendicular 
to the wall give a large contribution to 7// (-2) and a classical approximation is not adequate 
to describe the quantum interference nature of CP- violation. Secondly, the final result has 



the same infrared singularity responsible for the failure of the perturbative expansion 



for the values of the parameters such that rrih = (here the singular behaviour is smeared 
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out by the presence of the additional " mass " term ~ r/2). However, the theory remain 
perturbative as far as ah = {g^ l2-n){T'^ /m\) < 1, which gives nih ^ 0.3T. An estimate of the 
Higgs damping rate in the SM was obtained in Ref. [1^ in the low momentum limit and can 



be used here only to give very crude estimate of the Higgs coherent time, r ~ (10 — 10^)/T. 
With such values, our derivative expansion is barely justified since the wall thickness can 
certainly span the same range (10 — 10^)/T. We will discuss the validity of our approximation 
in the next section. 

The next step amounts to solving the set of coupled differential equations describing 
the effects of diffusion, particle number changing reactions and CP-violating source terms. 
We suppose that, among the supersymmetric particles, charginos, neutralinos and the right- 
handed stops as well as the Higgs degrees of freedom, are in equilibrium in the thermal bath. 
Under this hypothesis, strong sphalerons do not drive the asymmetry to zero [|l^. 

Closely following the approach taken in Ref. we can estimate the final baryon asym- 
metry generated by the Higgs sector to be 

(nB\ .ADT^^ 

— =-g{ki) — ^ , (6) 

\ s J H v^s 

where s = 27r^5f=KsT^/45 is the entropy density {g^,s being the effective number of relativistic 
degrees of freedom); g{ki) is a numerical coefficient depending upon the light degrees of 
freedom present in the thermal bath, D is the effective diffusion constant, T^g = Gna^rT is 
the weak sphaleron rate (/t ~ 1) ||19[0 and 



where 

2D ^ ^ 

The quantity T is the effective decay constant and 7_h-(z) = VwdzJfj{z)f{ki) is now defined 
in the bubble wall frame, /(fcj) being a coefficient depending on the number of degrees of 
freedom present in the thermal bath and related to the definition of the effective source [0. 
Since for relatively low values of the pseudoscalar mass (m^ > mz) the variation of the 
ratio of the vacuum expectation values of the Higgs fields along the bubble wall is small 



The correct value of k is at present the subject of debate, see, for instance, Ref. |2C] 
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@J2T|, we can take ^ ~ /5. For typical values tan/3 ~ 2, T ~ 5 x IQ-^T, ruh ~ 0.4 T and 
^vl{T) + vl{T)/T ~ 1.2, we find 



where the dependence upon v^, is very weak. 

The Higg scalar sector may account for the baryon asymmetry for large values of the 
phase (f)^. Such large values may be tolerated in view of the experimental limits on the 
neutron electric dipole moment of the neutron if the squarks of the first and second generation 
have masses of the order of a few TeV. 

3. One should not claim victory too soon, though. Let us look back at the approxima- 
tions we have made and discuss to which extent our estimate might be altered by performing 
a more refined analysis: 

a) Our result is rigously valid only for L^F > 1. Given our poor knowledge of both 
parameters, this condition does not seem unreasonable. Moreover, for L^F <^ 1 one may 
expect that the dependence of the final result upon F disappears so that J(/5, rrih, F) will 
saturate at the value F ~ L~^. Our results seem to confirm this expectation. 

h ) We have seen that quasiparticles with long wavelengths give the largest contribution 
to the source. This means that the classical approximation is not adequate to describe the 
quantum interference nature of CP- violation and a quantum approach must be adopted to 
compute the CP-violating source. This was done in the present paper. However, in the 
second stage of the computation, we have made use of classical Boltzmann equations. For 
low momentum particles, the validity of the classical Boltzmann equation starts to break 
down. It is indisputable that the ultimate answer can be provided only by a complete 
nonequilibrium quantum field theory approach. Kinetic theory and classical Boltzmann 
equations have been used to describe the dynamics of particles treated as classical with 
a defined position, energy and momentum. This requires that, in particular, the mean 
free path must be large compared to the Compton wavelength of the underlying particle 
in order for the classical picture to be valid, which is not guaranteed for particles with a 
small momentum perpendicular to the wall. Distribution functions obeying the quantum 
Boltzmann equations are the only correct functions to describe particles in an interacting, 
many-particle environment. 

Solving the quantum Boltzmann equations represents an Herculean task. However, recent 
investigations have revealed that the quantum Wigner distributions posses strong memory 




(9) 
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effects and that their relaxation time is typically longer than the one obtained in the classical 



limit 1 22]. This slowdown of the relaxation processes may keep the system out of equilibrium 
for longer times and therefore enhance the final baryon asymmetry. 

c) When computing the Higgs source, we have been assuming the loopwise expansion 
being reliable, in other words we have computed only the lowest order loop diagrams. How- 
ever, there is an infinite class of loop diagrams with more interaction vertices, the so-called " 
ladder " diagrams , which may contribute at the leading order even though full thermal 
corrections have been used in the evaluation of the lowest level diagrams. This is because 
there are additional infrared divergencies at finite temperature when the ladder- diagrams 
contain nearly massless Higgs modes and the infrared cut-off is provided by the mass rrih- 
The physical reason for this large correction has to be identified with the fact that the low- 
est loop diagram takes into account only the effects of particles which have undergo only 
a few collisions in the plasma, after which particles are not fully thermalized yet. Sum- 
ming up all the leading contributions is a nontrivial problem (in fact, it has been recently 
claim that the ladder graphs may be canceled by other types of terms in the loop expansion 
p^). A reasonable estimate of the relative size of of the contribution of the ladder-diagrams 
with respect to the lowest diagram computed in this paper is given by (loop suppression 
factor) X ((7^T^/m^) ln(T/m/i). Numerically this contribution may be sizeable and might 
change the final estimate of the baryon asymmetry by some factor 0{1). 
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